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We study ratios of partition functions in two types of sectors of 5(7(2), with fixed temporal center flux and 
with static fundamental charge. Both can be used as bona fide order parameters for the deconfinement transition. 



1. Introduction 

It is well known that the deconfinement tran- 
sition which occurs at a critical temperature Tc 
in Yang-Mills theory corresponds to the sponta- 
neous breaking of center symmetry. It is natural 
then to expect that the center degrees of free- 
dom of the gauge field are somehow responsible 
for confinement and its disappearance. Indeed, 
numerical studies tend to show that the center 
vortices, extracted from the gauge field after par- 
tial gauge-fixing, confine linearly at T ~ with 
the Y-M string tension at least in the contin- 
uum limit [Q. Here, we show quantitatively, in 
a gauge-invariant way, that the free energy of a 
center vortex is an order parameter for the decon- 
finement transition, as anticipated in . It jumps 
from to +00 as the temperature is raised. 

Alternatively, one should be able to use the free 
energy Fq of a static fundamental charge, jump- 
ing from -1-00 to at Tc. The Polyakov loop P is 
commonly used for this purpose in lattice studies. 
If (TrP) ~ e^T^i^ the symmetric (broken) phase 
gives for Fq an infinite (finite) value. However, 
the periodic boundary conditions (b.c.) within 
which (TrP) is measured are incompatible with 
the presence of a single charge. And, like any 
Wilson loop, (TrP) is subject to UV-divergent 
perimeter terms, such that (TrP) = at all tem- 
peratures as the lattice spacing a ^ 0. Here, fol- 
lowing Q, we measure the gauge- invariant, UV- 
regular free energy of a static fundamental charge, 
and show that it has the expected behaviour, dual 
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to that of a center vortex. Moreover, the string 
tension below Tc and the dual (vortex) string ten- 
sion above Tc Q are quantitatively related by a 
universal ratio, for more details, see 

2. 't Hooft loop, center vortices, and t.b.c. 

The 't Hooft loop W, defined on the dual lat- 
tice, multiplies every Wilson loop linked to it by 
a Zpf center element. In the case of an x Ry 
loop in SU{2), it is constructed by flipping the 
sign of the coupling /3 ^ — /3 of the Rx x Ry 
(z, i)-plaquettes dual to S, the rectangle bounded 
by W. [To compute the 't Hooft loop, we ac- 
tually flip f3 one plaquette at a time as in |^.] 
Equivalently, as shown explicitly by a change of 
integration variables, one could choose for E any 
other surface bounded by W. In any case, a set 
of plaquettes will carry the extra flux of a 
center vortex, whose precise location is gauge- 
dependent. The gauge- invariant 't Hooft loop 
W therefore generates a fluctuating center vor- 
tex sheet bounded by W. 

When an L X L spatial 't Hooft loop is created in 
an L'^ X Nt lattice, every L x Nt temporal plane 
dual to it must carry a flux. This is equivalent 
to enforcing twisted boundary conditions (t.b.c.) 
in this temporal plane. Calling Zq the periodic 
b.c. partition function, we measure the free en- 
ergy cost e"^*" = Zk{i)/Zo of t.b.c. in i = 1, 2 and 
3 temporal planes as a function of temperature T. 
Our results for Zk{l)/Zo are shown in Fig. 1 for 
lattices of size Nt = A,L = 6. .20. As expected 
from Finite-Size Scaling (FSS), they are well de- 
scribed around Tc by a single function of L/S^, 
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where the correlation length ^ diverges near Tc as 
with t = ^ — 1 the reduced temperature 
and V Ri 0.63 for the 3d Ising model universal- 
ity class. This is shown in Fig. 2, where the FSS 
variable is a; = LTc\t\'^. The large- cc asymptotic 
behaviour shows an area law exp{—aox'^) in the 
deconfined phase, and screening exp{exp{—bx)) 
in the confined phase. In the thermodynamic 
limit, Zk/Zo is a step function and can be used 
as a UV-regular order parameter as advertised. 
The interpretation of this behaviour is simple. 
Twisted b.c, or equivalently large center vortices, 
can be imposed at exponentially small cost when 
the system size becomes large compared to the 
confinement correlation length at T '--^ 0. At fi- 
nite temperature, the spread of the center vortex 
in a temporal plane is now limited by the com- 
pact dimension 1/T. Above Tc, the Polyakov loop 
acquires a non-zero expectation value. Tempo- 
rally twisted b.c, which enforce antiperiodicity 
on the Polyakov loop, then create an interface of 
minimal area separating positive and negative 
Polyakov loops. Such interfaces and the corre- 
sponding center vortices are exponentially sup- 
pressed by the dual-string or interface tension. 
T.b.c. in 2 and 3 temporal planes modify the min- 
imal area of the interface to and L^-^S re- 
spectively. We observe that the dual string ten- 
sion CTo is multiplied by precisely V2 and VS. 

3. Confinement/screening of fundamental 
charge: electric flux 

As seen above, temporally (say, {z,t)) twisted 
b.c. enforce antiperiodicity on the Polyakov loop: 
TrP^x)P{x + Lz) = -1, instead of +1 with pe- 
riodic b.c. Consider then the expectation value 
of this binary observable, in an enlarged parti- 
tion function Zeni which sums over all possible 
temporal b.c, periodic or twisted in each plane. 
Zcrd = Zo + 3Zfc(l) + 3Zk{2) + Zfc(3). Contrary 
to Zq, Z^ni allows for the consistent introduction 
of a static fundamental charge, which is neutral- 
ized by its image in the twisted direction (say z). 
Keeping track of the (z, i)-flux for each b.c. 

{TrPHx)P{x + Lz))=Z,il) (1) 
= J-(Zo + ZkH) - Zk{2) - Zfe(3)) 
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Figure 1. Partition function of one temporal twist 
as a function of temperature for various volumes. 

Similar expressions (in fact, Z2 Fourier trans- 
forms of the Zfe's) follow for 2 and 3 orthogo- 
nal 'electric' fluxes Fig. 3 shows ^^(1). Its 
behaviour is dual to that of Zfe(l). It also is 
a function of the FSS variable x, and it shows 
screening exp{exp{—bx'^)) above Tc and confine- 
ment exp{—aQx) below Tc, see Fig. 4. Moreover, 
because Z^ and Zk are simply related by Fourier 
transform, one obtains that b — and b = a^. 

Universality states that the FSS functions mea- 
sured here for SU{2) are the same as for the 3d 
Ising model after a single rescaling of x. This 
allows a precise determination of the large-x co- 
efficients (To, ctq by matching only the small-x part 
of the FSS functions between the two models 1^. 
Another universal ratio is given by 

a,/al = a{T+)/ {^^jpy « 0.41 . (2) 

Our data is consistent with this prediction, which 
relates the two sides of the phase transition. 

4. Summary 

Let us give an intuitive description of the de- 
confinement transition in Yang-Mills theory: 

• At low temperature, center vortices can spread 
to lower their free energy. Their proliferation dis- 
orders the Wilson loop and makes the average 
Polyakov loop traceless, reflecting confinement. 

• As the temperature is raised, the squeezing 
of center vortices in the compact 1/T dimen- 
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Figure 2. Free energy of one temporal twist as a 
function of the FSS variable x ^ for T ^ T^. 



Figure 4. Free energy of one electric flux as a 
function of the FSS variable x ^ for T ^ T^. 
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Figure 3. Partition function of one electric flux. 

sion prevents them from arbitrarily lowering their 
free energy. Eventually, they cannot disorder the 
Polyakov loop anymore. Macroscopic regions of 
Polyakov loops of a definite center sector appear, 
separated by interfaces whose tension exponen- 
tially suppresses center vortices. Confinement has 
disappeared. This squeezing does not affect cen- 
ter flux in spatial planes, so that the spatial string 
tension persists at all temperatures. 

Because the free energy Fq of a fundamental 
charge results from the cancellation of partition 
functions with different b.c, — e~~^« is expo- 
nentially small in the confined phase, and ~ 1 in 
the deconfincd phase. We have shown that the ra- 
tios of partition functions Zk/Zg or Z^/Zq, which 
are well-defined continuum quantities, can both 
be used as order parameters for the transition. 



If dynamical fundamental charges (quarks) are 
introduced in the system, the picture becomes 
murky, reflecting the change from a phase transi- 
tion to a crossover. To start with, the 't Hooft 
loop now depends on the particular surface S 
along which center flux is introduced. Its rele- 
vance in full QCD remains to be studied. 

We thank M. Garcia Perez, M. Hasenbusch and 
M. Pepe for discussions. 
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